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Abstract 

^ . We compute the effective hamiltonian for non-leptonic | AF\ = 1 decays in the standard 

model including next-to-next-to-leading order QCD corrections. In particular, we present 
the complete three-loop anomalous dimension matrix describing the mixing of current- 
current and QCD penguin operators. The calculation is performed in an operator basis 
which allows to consistently use fully anticommuting 75 in dimensional regularization at 
an arbitrary number of loops. The renormalization scheme dependences and their can- 
cellation in physical quantities is discussed in detail. Furthermore, we demonstrate how 
our results are transformed to a different basis of effective operators which is frequently 
adopted in phenomenological applications. We give all necessary two-loop constant terms 
which allow to obtain the three-loop anomalous dimensions and the corresponding initial 
conditions of the two-loop Wilson coefficients in the latter scheme. Finally, we solve the 
renormalization group equation and give the analytic expressions for the low-energy Wilson 
coefficients relevant for non-leptonic B meson decays beyond next-to-leading order in both 
renormalization schemes. 



1 Introduction 



Perturbative QCD effects have an important impact on the structure of the effective hamil- 
tonian for non-leptonic \ AF\ = 1 processes with F = S, C or B, which describes the weak 
decay of the corresponding mesons and hadrons. Most notably, they can lead to a sizable 
enhancement of the AI =1/2 transitions, of the CP- violating ratio e'/e, and of the QCD 
penguin contributions to rare and radiative B decays within the Standard Model (SM) [1] 
and some of its innumerous extensions [2-7]. 

In all cases, these short-distance QCD effects can be systematically calculated using 
an effective field theory framework, which allows to resum large QCD logarithms of the 
form L = \nfi/M w by solving the Renormalization Group Equation (RGE) that governs the 
scale dependence of the Wilson coefficient functions of the relevant \AF\ = 1 local operators 
built out of the light and massless SM fields. After the pioneering Leading Order (LO) 
calculation of the 0(o;"L n ) contributions [8], the resummation of the 0{a 1 ^L n ~ 1 ) logarithms 
has been completed more than ten years ago and subsequently confirmed by several groups. 
The main components of the perturbative Next-to-Leading Order (NLO) calculation are 

i) the one-loop 0(a s ) corrections to the relevant Wilson coefficient functions [9-11] and 

ii) the two-loop 0(a 2 ) Anomalous Dimension Matrix (ADM) describing the mixing of the 
associated physical operators [10-15]. 

To improve on the present NLO calculation, one needs to include one more order in 
the strong coupling expansion, aiming at a resummation of all the 0(a™L n ~ 2 ) logarith- 
mic enhanced corrections. The completion of this Next-to-Next-to-Leading Order (NNLO) 
computation constitutes the core of this work. Since the two-loop 0{a 2 s ) matching cor- 
rections to the relevant Wilson coefficients are already known from [16] the only missing 
ingredient to perform this task is the knowledge of the three- loop 0(al) ADM describing 
the mixing of the current- current and QCD penguin operators. In this paper we will close 
this gap by employing standard techniques [15, 17, 18] to carry out a direct calculation of 
the required ADM adopting the renormalization scheme introduced in [11]. Since in this 
scheme the QCD penguin operators are defined in such a way that traces with 75 do not 
occur to all orders in perturbation theory, we are allowed to consistently use dimensional 
regularization with a naive anticommuting 75. This feature is very welcome, as it makes 
the actual three-loop calculation completely automatic and rather straightforward. 

The NNLO ADM we have computed can be used in analyses of new physics models as 
well, provided they do not introduce new operators with respect to the SM. This applies, for 
example, to the case of the two Higgs doublet models [2], to some supersymmetric scenarios 
with minimal flavor violation [3], and to specific models of universal extra dimensions [4]. 
On the other hand, in left-right-symmetric models [5] and in the general supersymmetric 
SM [6], additional operators with different chirality structures arise [19]. In many cases 
one can exploit the chiral invariance of QCD and use the same ADM, but in general an 
extended basis is required. 
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Another strong motivation to write the article at hand was that the three-loop ADM 
computed here is part and parcel of the complete NNLO analysis of rare semi-leptonic 
B — > X s £ + £~ decays presented recently by us in collaboration with Christoph Bobeth and 
Paolo Gambino [20]. Furthermore, it constitutes a integral part of the NNLO calculation 
of radiative B — > X s ^ decays, admittedly a very ambitious enterprise, which nevertheless 
has already aroused the interest of some theorists [21]. Whereas our previous work dealt 
exclusively with the phenomenological application of our result, we will now — in the 
spirit of [9-11] - focus on the more formal aspects of the renormalization of effective 
field theories, such as the issue of scheme dependences in general, their cancellation in 
physical observables, and the general transformation properties of the ADM and the Wilson 
coefficients under a change of scheme. In this respect we will extend the existing NLO 
results [9,10] to the next order, paying special attention to the conceptual features related 
to the renormalization of the strong coupling constant. 

While much of the discussion in this paper is therefore rather technical, our general 
results have important practical applications. This will be showcased by means of a couple 
of examples. In particular, we will devote a sizable part of the present article to derive 
the explicit NNLO relation between our and a different renormalization scheme that is 
commonly used in the literature on weak decays [1]. In the latter scheme, which we shall call 
"traditional" scheme from now on, the applied form of the effective hamiltonian introduces 
unwanted traces containing 75 by definition. These traces turn out to be harmless at the 
LO, but involve a lot of technical difficulties related to the use of Fierz symmetry arguments 
in n = 4 — 2e dimensions at the NLO [10, 14]. Applying the same scheme in a direct 
calculation of the ADM at the NNLO or even beyond would thus be extremely tedious. 
We will not attempt such a direct computation here, but avail the derived NNLO relation 
between our and the "traditional" scheme to find the NNLO ADM and the corresponding 
matching conditions in the latter scheme on detours. As another exercise we solve the RGE 
and give the analytic expressions for the AB = —AS = 1 low-energy Wilson coefficients 
beyond NLO in both renormalization schemes. We are aware of the fact that some of the 
formulas presented below are rather long. Nevertheless we believe that at least some of 
them should be useful to the reader interested mainly in the application of the presented 
formalism to weak decays rather than in the conceptual subtleties, which obviously address 
more technical minded colleagues. 

The main part of this paper is organized as follows: in Section 2 we present the general 
structure of the effective hamiltonian for non-leptonic \AF\ = 1 decays at the NNLO 
level. Section 3 is devoted to the simplest application of the general formalism, namely 
the AB = —AS = 1 decays. We recall the relevant effective hamiltonian and list all the 
dimension-five and six operators that will be needed in the calculation of the three-loop 
ADM. In Section 4 we collect the results for the initial conditions of the relevant Wilson 
coefficients through NNLO obtained in [11, 16]. After a brief description of the actual 
three- loop calculation the final result for the ADM is given in Section 5. In Section 6 we 
solve the RGE to find the explicit NNLO expressions for the low-energy Wilson coefficients 



3 



relevant for non-leptonic B meson decays. In Section 7 we elaborate on the question of 
scheme dependence related to the renormalization of the effective operators as well as the 
strong coupling constant. Section 8 starts out with a general discussion of the non-trivial 
nature of a change of the basis of physical operators in the framework of dimensional 
regularization, followed by a demonstration of how the NNLO results for the ADM and 
the matching conditions are transformed to the "traditional" basis of effective operators. 
Finally, in Section 9 we summarize the main results of this work. 

Some technical details as well as additional material has been relegated to the appen- 
dices: in Appendix A.l we derive the explicit form of the matrix kernels that are needed 
to find the evolution matrices through NNLO, while Appendices A. 2 and A. 3 contain all 
ingredients that are necessary to transform our results to the "traditional" set of oper- 
ators. The NNLO analytic formulas for the low-energy Wilson coefficients relevant for 
non-leptonic B meson decays in the latter scheme will be given in Appendix A. 4, which 
concludes our paper. 



2 General Structure 

The effective hamiltonian for non-leptonic |AF] = 1 decays has the following generic 
structure [1] 

n cS = --^-v CKM Q T c(fi). (i) 

Here G F denotes the Fermi constant and Q T is a row vector containing the relevant local 
operators Qi. Explicit expressions will be given in Section 3. C(fi) is a column vector con- 
taining the Wilson coefficients Cj(/i) that together with the Cabibbo-Kobayashi-Maskawa 
(CKM) factor [22] Vckm describe the strength with which a given operator enters the 
hamiltonian, and [i is the renormalization scale. The decay amplitude for a decay of a 
meson M into a final state F is simply given by (F\H c g\M). 

The Wilson coefficient functions evolve from the initial scale /i down to /i, which in 
practical applications is much lower than /io, according to their RGE. Using dimensional 
regularization with n = 4 — 2e and considering only mass independent renormalization 
schemes it is given by 

^CQm) = i T (g)C^) , (2) 

where j(g) is the ADM corresponding to Q. Neglecting all electromagnetic effects, the 
general solution of this equation reads 

C(fi) = U(^fio)C(fx ), (3) 
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with 

U(vL,iM>) = T g exp dg'^-, (4) 

°° / a 2 °° / a 2 

T<", and = ^ A, (5) 

j=0 V 7 i=0 v 7 

Here C(/io) are the initial conditions of the evolution and T g denotes ordering of the coupling 
constants g(fx) in such a way that their value increases from right to left. (3(g) is the QCD 
(3 function. 

Keeping the first three terms in the expansions of 7(g) and (3(g) as given in Eq. (5), we 
find for the evolution matrix U (fi, /i ) in the NNLO approximation 

U(fi,fi ) =K( f i)U (0 \ f i, f i )K- 1 ( f x ), (6) 

where 



and 



k(^=i + ^jv + (^y.p\ 

k-\vo) = 1 - - (^) 2 (y (2) - 



(7) 



(8) 



denotes the LO evolution matrix, which depends on the matrix V and the so-called magic 
numbers a, that are obtained via diagonalizing ■yW 1 ': 

(V-^ T V) = 2/W^. (9) 

V / ij 

In order to give the explicit expressions for the matrices and we define 

jd) = VS^V' 1 , and 6 {i) = V-^^V , (10) 
for i = 1, 2. The entries of the matrix kernels and are given by 



where the first line recalls the classical NLO result [9], and the second one represents the 
corresponding NNLO expression, in agreement with [23]. The explicit derivation of S^ 1 ) 
and is presented in Appendix A.l. 

Let us now recall how the initial conditions of the Wilson coefficients are obtained. 
The amplitude for a given non-leptonic quark decay is calculated perturbatively in the full 
theory including all possible diagrams such as W^-boson exchange, QCD penguin and box 
diagrams as well as gluon corrections to all these building blocks. The result up to the 
NNLO is given schematically by 

Auii = (Q) (0)T (a^ + + , (12) 

where {Q) ^ denotes the tree-level matrix elements of Q. 

A second step involves the calculation of the decay amplitude in the QCD effective 
theory. It generally requires the computation of the operator insertions into current- current 
and QCD penguin diagrams of the effective theory together with gluon corrections to these 
insertions. Including QCD corrections up to the NNLO one finds 

A cS = (Q)^ T (l + ^Mf^ T + (^)V )T ) CM , (13) 

where the matrices f ^ and f^> codify the one- and two-loop matrix elements of Q, respec- 
tively. 

The matching procedure between full and effective theory establishes the initial condi- 
tions C(hq) for the Wilson coefficients. Equating Ami and A e s in Eqs. (12) and (13) at a 
scale hq translates into the following identity [24] 



+ ^ gfOfjO y ^(2) _ -(1)T _f(l)T£(p) _^(2)T^(0)j _ 

Combining Eqs. (2), (6), (7) and (14), we finally obtain 

cm = k(/^ o Wo) ( i (o) + - ^ (i) ^ (o) 

\ 4"7T 



(14) 



+ 

where 



47T 



(15) 



R W = f^ T + jW , and m = f^ T + + r« T J« , (16) 

are certain combinations of r^ T , f^ T , and which will play a special role in 
Section 7 where we will discuss the issue of renormalization scheme dependences in detail. 
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3 Effective Hamiltonian for AB = —AS = 1 Decays 



The simplest application of the general formalism outlined in the previous section is the 
case of non-leptonic B meson decays governed by the b — > s transition. For definiteness 
we will therefore give explicit formulas for the AB = —AS = 1 decays only. However, it is 
straightforward to transform them to the other |AF| = 1 cases. Neglecting contributions 
proportional to the small CKM factor V* s V u b which are irrelevant here, the corresponding 
effective off-shell hamiltonian is given by 

H cS = -^£-V*V tb (Q T C(fi) + N T C N (fi) + B T C B {ti) + E T C E ^)) . (17) 

In addition to the gauge-invariant operators Q, non-physical operators arise as countert- 
erms in the renormalization of higher loop One-Particle- Irreducible (1PI) off-shell Green's 

— * 

functions with insertions of the operators Q. These non-physical operators can in general 
be divided into three different classes [13,25-27]: i) operators N that vanish by use of the 
QCD Equations Of Motion (EOM), ii) non-physical counterterms B that can be written 
as a Becchi-Rouet-Stora-Tyutin (BRST) variation [28] of some other operators — so-called 
BRST-exact operators — and Hi) evanescent operators E that vanish algebraically in n = 4 
dimensions. 

— * 

The set of physical operators Q consists of six dimension-six operators, which can be 
chosen as [11,29] 



Qi = 




Q2 = 


(s L 'Yv 1 CL)(c L 'f 1 b L ) , 


Q 3 = 




Qa = 


(s L ^T a b L )J2 (Ql^q), 


Q 5 = 






{sLl m T% L ) (qi m T a q) , 



(18) 



and one dimension-five operator 

Qs = -m b (s L <j^T a b R )G a ^ 2 . (19) 

Here we have used the definitions 7 /n ... /lm = 7 M1 • • -7 Mm , 7 M1 "' Mm = 7 Ml • • -7^™ and a^ 1 ^ 2 = 
i [7 Ml , 7 M2 ]/2, and the sum over q extends over all light quark flavors, g is the strong 
coupling constant, qi and qn are the chiral quark fields, G^ 1At2 is the gluonic field strength 
tensor, and T a are the generators of SU(3) C , normalized so that Ti (T a T b ) = 5 ab /2. 
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Figure 1: The one-loop 1PI diagrams which mix Q 2 into N\ 



(i) 



The physical operators given in Eqs. (18) and (19) include the current-current operators 
Qi and Q 2 , the QCD penguin operators Q^-Qs and the chromomagnetic moment operator 
Qs- Notice that we have defined Qi~Q& in such a way that problems connected with the 
treatment of 75 in n = 4 — 2e dimensions do not arise [11]. Consequently, we are allowed 
to consistently use a fully anticommuting 75 in dimensional regularization throughout the 
calculation. 

As far as the EOM- vanishing operators are concerned, the specific structure of only one 
of them [11] 

1 



N\ L > = -s Ll ^T a b L D^Gl lil2 + Q 4 , 



is relevant in finding the one- and two-loop mixing of the four-quark operators Q\-Q%. The 
corresponding divergent one-loop 1PI diagrams are shown in Figure 1. 

In order to remove the ultraviolet (UV) divergences related to the two-loop subdiagrams 
with insertions of Q\-Q§ depicted in Figure 2, another ten EOM-vanishing operators need 
to be considered [15] 



1 



= 



iVf > = - 



m b s L ppb R , 



2 s L pppb L , 
[ 7s L p a^T a b L G a ^ 2 - G a ^ 2 7s L T a a^pb L 



+ Q 8 , 



N, 



(2) 



-m b s L 
9 



pG-Gp 



N, 



(2) 



s L \pGG- GGpj b L - im b s L GGb 



R 



(2) 



1 

9 



N- 



(2) 



s L (pp G + Gpp ) b L + im b s L Gpb R 



s L ( p G^G^ - G^G^p ) b L - im^G^G^bR 



Ml 
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Figure 2: Some of the two-loop 1PI diagrams which mix Q 2 into N[ 1] and n[ 2) -N$. 



ivf > = 1 - 



iv 9 (2) = 



s L (pD^ + G^p ) b L + im b s L G^D^b R 



N$ = d abc 



s L p Gpb L + im b s L p Gb R 
s L (p T a - T a p] b L - im b s L T a b R 



G^G^, 



(21) 



where = <9 M + igG^ and = — igG^ denotes the covariant derivative of the gauge 
group SU(3)c acting on the fields to the right and left, respectively. G a denotes the gluon 



field, and we have used the definitions G, 



G a t T a and d abc 



2Tr({T a ,T 6 }T c 



It is important to remark that the EOM-vanishing operators introduced in Eqs. (20) 
and (21) arise as counterterms independently of what kind of infrared (IR) regularization 
is adopted in the computation of the anomalous dimensions of Qi-Qe- However, if the 
regularization respects the underlying symmetry, and all the diagrams are calculated on- 
shell, non-physical operators have vanishing matrix elements [25-27,30]. In this case the 
EOM-vanishing operators given in Eqs. (20) and (21) play no role in the calculation of the 
mixing of physical operators. If the gauge symmetry is broken this is no longer the case, 
as diagrams with insertions of non-physical operators will generally have non-vanishing 
projection on the physical operators. As we will discuss in Section 5, our IR regularization 
implies a massive gluon propagator, and therefore non-physical counterterms play a crucial 
role at intermediate stages of the anomalous dimensions calculation. 

In contrast to the case of the two-loop mixing of the magnetic operators considered 
in [15,17], it is a priori not clear if BRST-exact operators do arise as counterterms of Q\- 
Qq. Since the BRST variation raises both ghost number and mass dimension by one unit, 
it is evident that any BRST-exact operator that potentially could mix with Q\-Qq has to 
be a BRST variation of a dimension-five operator containing a single anti-ghost field. The 
only possibility for the latter operator having the correct chirality structure is given in the 
Ri: gauge by 



= s 



1 

9 



jd^G^-gf^ {d^u h ) 



u 



{s Ll ^T a b L ), (22) 
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(a) (b) (c) 

Figure 3: (a) A typical example of a divergent two-loop 1PI diagram which potentially could 
introduce a mixing of Q2 into (b) A typical example of a counterterm contribution 

needed to renormalize the corresponding two-loop 1PI diagrams, (c) The contribution to 

(2) — 
the one-loop matrix element of B\ containing an effective b — > su a u a vertex, which has a 

non- vanishing on-shell projection on Q4 if a non-zero ghost mass is used in the calculation. 

The QCD penguin and box contributions to the matrix element that contain an effective 

b — > sg vertex are not shown. 



where s denotes the BRST operator, u a and u a are the ghost and anti-ghost fields, f abc 
are the totally antisymmetric structure constants of SU(3)c and £ is the covariant gauge- 
parameter. 

(2) 

Although there is no obvious reason why B\ should not appear as a counterterm of 

(2) 

Qi^Qqi it turns out that up to three loops B\ does not play a role in the mixing of physical 
operators considered in the paper at hand. The key observation thereby is that the overall 
contribution from the two-loop 1PI diagrams depicted in Figure 3 (a) is canceled by the 
corresponding counterterm contribution as shown in Figure 3 (6), so that the associated 
renormalization constant is exactly zero at 0(a 2 s ). Therefore B^ does not contribute to 
the mixing of Q\-Qq into Q4, although its one-loop 0(a s ) matrix element displayed in 
Figure 3 (c) does not vanish if it is computed using non-vanishing gluon and ghost masses 
to regulate IR divergences. 

In order to remove the divergences of all possible 1PI Green's functions with single 
insertion of Q\-Q% we have to introduce some evanescent operators E as well. At the 
one-loop level one encounters four evanescent operators, which can be chosen to be [11,29] 

E[ 1] = {s Llm T a C L ){c Ll ^T% L ) - I6Q1 , 

= (s Ll ^ s c L )(c Ll ^b L ) - 16Q 2 , 
Ei 1] = {s Llmm b L ) ^w^ g ) + 64Q 3 - 20g 5 , (23) 

E? = {s Llmm T a b L ) (qr im " 5 T a q) + 64Q 4 - 20Q 6 . 

At the two-loop level four more evanescent operators do arise, that can be defined as [11,29] 

E[ 2) = (s L 7 W M 2 M3M4/, 5 ^ a CL)(cL7 WM2M3M4M5 T a 6L) - 256Q! - 20E? , 
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= (sL7M 1 M 2 M3M4M 5 CL)(c L7 mM2M3M4M5 &L) - 256Q 2 - 204 1} , 
Ef ] = {s Lllxim „ m b L ) + 1280g 3 - 336g 5 , 

Ef = (s L ^ imm T a b L ) J2 q (qr imm T a q) + 1280Q 4 - 336Q 6 . (24) 

Finally, at the three-loop level another four evanescent operators are needed. We define 
them in the following way: 

E[ 3) = (sL7 W M 2W /,4/ 15W ^^CL)(cL7 WAt2 ' 13At4 ' 15At6 ' 17 T a 6 L ) - 4096Q! - 336E; i} , 

E { 2 ] = (s L7w M 2W M4^M6/ 17 CL)(cL7 WAt2 ' 13At4 ' 15At6 ' 17 ^) - 4096Q 2 - 3364 1} , 

£f = {sLl mmm M J2 q (qY imm ^ 9 q) + 21504Q 3 - 544og 5 , (25) 

Ef ] = (s Llmim6m T a b L ) J2 q (qr mm, * m T a q) + 21504Q 4 - 5440Q 6 . 

(3) (3) 

Needless to say, the above choice of evanescent operators E\ -E\ is not unique, in the 
sense that their particular structure can be changed quite a lot without affecting the three- 
loop ADM of the four-quark operators Q\-Q&. For instance, adding any multiple of e times 
any physical operator to them leaves the ADM unchanged up to 0{a^). This is in contrast 
to what happens if such a redefinition is applied to the one- and two-loop evanescent 
operators given in Eqs. (23) and (24). 



4 Initial Conditions of the Wilson Coefficients 

Let us now turn to the initial conditions C(/xo) of the Wilson coefficients. Their values are 
found by matching the full to the effective theory amplitudes perturbatively in a s . The 
NLO and NNLO approximation requires the calculation of one- and two-loop diagrams both 
in the SM and the low-energy effective theory. Some of the SM two-loop 1PI diagrams one 
has to consider in order to find the 0(a 2 s ) corrections to C(yU ) are displayed in Figure 4. 
Restricting ourselves to the physical on-shell operators Qi-Qe and setting /x = M w , one 
obtains using dimensional regularization with a naive anticommuting 75 [11,16]: 

nti*\ , c a»(A/»0 , f a,(M w ) \ 2 [7987 17 2 ~ \ 

C 2( M^ 1+ (^») 2 (if + V 
C3(M„)=(^) 2 g l( x ( ), 
Ci (M w ) = ^E ( Xt ) + (^>Y £,(,,) , 
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b b W 



Figure 4: Some of the SM two-loop 1PI diagrams one has to calculate in order to find the 
Wilson coefficients of the four-quark operators Q\-Q% at 0(a 2 ) 



C e (M w ) = g + \e«(z<) - , (26) 

where x t = mf/M^. The one-loop Inami-Lim [31] function E (x t ) characterizing the 
effective off-shell vertex involving a gluon reads 

~ 8 - A2x t + 35a; 2 - 7x 3 t 4 - 16x t + 9x 2 

En(Xt) = : rr; : ri mx t ■ 27 

V ; 12(x t -l) 3 6(x t -l) 4 V ; 

The one-loop function T (x t ) originates from diagrams like the first one shown in Figure 4. 
Subtracting the corresponding terms in the gluon propagator in the so-called momentum 
space subtraction scheme at q 2 = 0, which guarantees that a s has the same numerical value 
on the full and effective side at the matching scale through NNLO, one finds [16] 

~. . 112 /20 \ 

To{x t ) = — + 32x t + — + lox t hvx t 

(i w (28) 

- (8 + 16x t ) V^^l Cl 2 \2 arcsin [iq= J J , 

with Ch(x) = lm[Li 2 (e lx )} and Li 2 (x) = — J* dt ln(l — t)/t. The remaining two-loop 
functions E\{x t ) and G\(x t ) take the following form [16] 

1120 - 12044x t - 5121^ - 5068a; 3 + 7289a;? 



El(Xt) - 648(s t -l)' 

380 - 7324x t + 17702a; 2 + 2002a;, 3 - 5981a;? + 133x5 , 

lnx t 



324(x t - l) 5 

112 - 530x t - 3479a; 2 + 2783a: 3 - 1129a;? + 515x5 , 2 

—, 1 hi x t 

108(x t - l) 5 

40 - 190a; t - 81a;, 2 - 614a:? + 515a;? . , 10 , 
54(x^l7 Wl + , 
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~ _ 554 - 2h23x t + 2919x t 2 - QQ2x 3 t 243(x t - l) 3 

88 - U2x t - 357x1 + 100x 3 + 35xf _ 20 - AOx t + hxj 2 



11(1, -l) 4 27(„ -i) 2 

40 - 160,, - 3te ? + 100,° - 10,? _ _ 20 

27(:r t -l) 4 2V ^ 81 

5 Anomalous Dimension Matrix 



Before presenting our results for the anomalous dimensions describing the mixing of the 
four-quark operators Q\-Qq up to 0(a^) let us recall some definitions that will turn out 
to be useful in the rest of the paper. 

Upon renormalization the bare Wilson coefficients Cb{^) of Eq. (1) transform as 

C B (fi) = Z T C(fx). (30) 
In terms of the renormalization constant matrix Z the ADM of Eq. (2) is then given by 

7( 9 ) = ^-'. (3D 

The renormalization constants Zij of the operator Qj can be expanded in powers of g in 
the following way 



k=l ^ ' 1=0 



(32) 



Following the standard MS scheme prescription, Z^ is given by pure l/e l poles, except 
when % corresponds to an evanescent operator and j does not. In the latter case, the 
renormalization constant is finite, to make sure that the matrix elements of the evanescent 
operators vanish in n = 4 dimensions [13,27]. 

In a mass independent renormalization scheme the only /z-dependence of Zy resides in 
the coupling constant. In consequence, we might rewrite Eq. (31) as 

la = 0(e, 9) Z ik^j-Z-f , (33) 

where (3(e,g) is related to the usual QCD (5 function via 

P(e,g) = -eg + P(g). (34) 
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Figure 5: Some of the three-loop 1PI diagrams we had to calculate in order to find the 
mixing among the four-quark operators Qi~Qe at 0(af). 



The finite parts of Eq. (33) in the limit of e going to zero give the anomalous dimensions. 
Inserting the expansions of 7(g) and j3{g) in powers of g, as given in Eq. (5), one im- 
mediately finds [15, 18] for the anomalous dimensions governing the evolution of physical 
operators up to third order in the strong coupling parameter: 

7 (°) = 2Z^ , 

7 M = AZ^ - 2Z^Z^ , (35) 
^(2) = 6Z (3>1) - 4Z {2 ' 1] Z (lfi) - 2Z (1 ' 1) Z {2 ' 0) . 

The matrices Z (lfi \ Z^\ Z (2 ^ and Z^ are found by calculating various one- and 
two-loop diagrams with a single insertion of Qi~Qe, E^-E^ and E[ 2 ^-E^\ whereas 
the matrix Z^ 3 ' 1 ^ requires the computation of three-loop diagrams with insertions of Q\- 
Qq as shown in Figure 5. The pole and finite parts of these one-, two- and three-loop 
diagrams are evaluated using the method we have described together with Paolo Gambino 
in detail in [15]: We perform the calculation off-shell in an arbitrary gauge which allows 
us to explicitly check the gauge-parameter independence of the mixing among physical 
operators. To distinguish between IR and UV divergences we follow [17,18] and introduce 
a common mass M for all fields, expanding all loop integrals in inverse powers of M. This 
makes the calculation of the UV divergences possible even at three loops, as M becomes 
the only relevant internal scale and three-loop tadpole integrals with a single non-zero 
mass are known [18,32]. On the other hand, this procedure requires to take into account 
insertions of the non-physical operators and N^-N$, as well as of appropriate 
counterterms of dimension-three and four, some of which explicitly break gauge invariance. 
A comprehensive discussion of the technical details of the renormalization of the effective 
theory and the actual calculation of the operator mixing is given in [15]. 

Having summarized the general formalism and our method, we will now present our 
results for an arbitrary number of quark flavors denoted by /. For completeness we start 
with the regularization- and renormalization-scheme independent matrix 7^, which is 
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given by 



7 



(0) 



V 



-4 


8 
3 





2 
9 








12 








4 
3 

















52 
3 





2 








40 

9 


160 ,4 f 
9 ^3 J 


4 
9 


5 
6 











256 
3 





20 








256 


544 , 40 f 
9 + 3 J 


40 


2 


9 


9 


3 



(36) 



While the matrix 7^ is renormalization-scheme independent, ■y^ and 7^ are not. In the 
MS scheme supplemented by the definition of evanescent operators given in Eqs. (23), (24) 
and (25) we obtain 



/ 



7 



and 



(i) 



35 \ 

1 R9 



/ 



7 



(2) 



145 , 16 r 
3 + 9 J 


-26+f 


f 1412 
J 243 


1369 
243 


134 
243 


35 
162 


-45+f / 


28 


416 


1280 


56 


35 


3 


81 


81 


81 


27 








4468 
81 


29129 52 f 
81 9 3 


400 
81 


3493 2 f 
108 9 J 








13678 , 368 f 
243 81 


79409 | 1334 r 
243 81 


509 8 p 
486 81 J 


13499 5 f 
648 27-' 








244480 160 <• 
81 9 J 


29648 2200 x 
81 9 J 


23116 1 16 f 
81 " r 9 J 


3886 1 148 r 
27 9 ■> 








77600 1264 f 
243 81 •< 


28808 , 164 f 
243 81 ■> 


20324 | 400 j- 
243 81 


21211 | 622 f 
162 " 1 " 27 •> 


1927 , 257 
2 ' 9 


f+ff+ 


(224+if/)C3 


475 , 362 
9 27 


r 40 r2 /" 896 

j~27-f ~l~r _ 





(37) 



5875184 
6561 



307 + 3|l / _^ / 2_ (1 34 4+160/)C 3 






269107 2288 f 1360 a 
13122 729 J 81 ^ 3 
69797 I 904 j, 2720 a 
2187 "•" 243 3 ' 27 ^ 3 
4203068 , 14012 f 608 a 
2187 ~<~ 243 3 27 ^ 3 
217892 . '-' 



1298 76 
3 



3-/-224C3 








2425817 , 30815 f 776 a 



13122 
1457549 



1 4374 J 81 
22067 f 2768 



8748 729 ■> 

18422762 , 888605 j- , 272 f 2 



C3 



317892 f , 472 r2 , f 27520 , 1360 f V 
6561 1 2187 81 3 ""v 81 " l ~ 9 3 I ' 3 

194951552 , 358672 j 2144 f 2 i 87040 a 
2187 ' 81 J 81 J ~ r 27 ^ 3 
162733912 2535466 j , 17920 j2 i / 174208 , 12160 a 
6561 2187 J " r 243 J ' V 81 ' 9 3 ) * 3 

r | 124 a 



2187 
70274587 



13122 
130500332 
2187 
13286236 



1 2916 
8860733 
17496 
2949616 



f 4010 y2 
3 729 J 



27 
f 39824 i 
I 27 1 

f 16592 

I 81 



160/) Cs 



343783 


, 392 


52488 


"'" 729 


37889 


28 


8748 


243 ■' 


674281 


1352 



--, J ,-/+ 3 §f/2 + (2 3 §M+640/)C 3 

1826023 f 159548 f 2 _ I 24832 , 9440 jAa 
3 729 J I 81 T 27 • / / S > 3 



^/+^C3 



4374 



69984 ^ 972 
366919 35 



4374 
31175 



243 
52 



52488 8748 ■> 81 
14732222 27428 



2187 81 
22191107 I 395783 
13122 t 4374 



/-f/ 2 -( 



27 

/"^C3 

3154 , 
81 " + 



11664 
9284531 2798 



__35r_ 11U A 
162 3 9 "» J 



136 
9 



272 f 2 13984 
81 ■< 27 



f 1720 r2_f 
J 243 3 \ 



81 



/)Cs 

±f>/)<3 



/"f/ 2 -( 1 f i +20/)C3 



11664 

3227801 105293 r 65 f 2 , I 200 _ 220 /A a 
11664 J 54 J 27 9 3 J^ A 



8748 
16521659 
2916 

32043361 I 3353393 
8748 ~ l ~ 5832 



, 4 /-^/ 2 -(^+200/)C 3 

/ _^3 / 2 + (9248_1^0 / ) C3 y 



(38) 



As far as the one- and two-loop mixing of the four-quark operators Q\-Q§, namely 7^ and 
7W are concerned, our results agree with those of [11]. Furthermore, they also agree with 
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the results obtained in [10, 14] after a transformation to the "traditional" operator basis. 
This will be shown in Section 8 by an explicit calculation. On the other hand, the three- 
loop mixing of Q\-Qq described by is entirely new and has never been given before. 
As it is characteristic for three-loop anomalous dimensions the entries of contain terms 
proportional to the Riemann zeta function £ 3 . 



6 Renormalization Group Evolution 

In this section we shall use the obtained ADM to find the explicit NNLO expressions for 
the Wilson coefficients 



cM) = c?\^) + u -^'' ll] <-^ inA '"' ! ^ 



4:71 



ill 



(39) 



with i = 1-6, at the low-energy scale fib = 0(mb), which is appropriate for studying non- 
leptonic B meson decays. Using the general solution of the RGE given in Eq. (6), we arrive 
at 



a 



s(o) 



6 



6 



(2) , (2) 2 



Eo(x t ) 



where r\ = a s (M w ) / a s (fib) and 



( 23 ~M °- 4086 -0-4230 -0.8994 0.1456 ) , 



/ 1 -1 



2 
3 

_2_ 

63 
J_ 

21 



1 

3 

__\_ 

27 
1 
9 

1 1 

126 108 

J_ __L 

84 36 







-0.0659 0.0595 -0.0218 0.0335 

0.0237 -0.0173 -0.1336 -0.0316 

0.0094 -0.0100 0.0010 -0.0017 

0.0108 0.0163 0.0103 0.0023 / 



(40) 



(41) 



(42) 
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5.9606 
1.9737 
-0.5409 
2.2203 
0.0400 
\ -0.2614 



/ 



2.0394 
1.3596 
0.0647 
0.0971 
-0.0162 
-0.0243 



1.0951 
-1.3650 
1.6332 
2.0265 
-0.1861 
-0.1918 

5.9049 
-1.9683 

0.2187 
-0.6561 
-0.0547 

0.1640 






1.6406 
-4.1830 
-0.1669 

0.4197 




-0.4268 
0.1534 
0.0606 
0.0700 






-1.6702 
-0.7135 
0.1887 
0.0295 




-0.5165 
0.1500 
0.0865 

-0.1412 



























-0.1933 








0.0695 








0.0274 








0.0317 






0.1579 
-0.0459 
-0.0264 

0.0432 






0.1428 
0.8752 
-0.0064 
-0.0675 






-0.2576 
-1.8215 
0.0201 
0.1474 




0.2832 
1.7355 
-0.0128 
-0.1339 

\ 



\ 


-0.2250 

0.7996 

0.0304 
-0.0640 / 

\ 


-0.2034 

0.1916 

0.0103 
-0.0140 / 






-0.1074 

0.1012 

0.0055 
-0.0074 j 



( 56.4723 
14.7825 
1.9905 
8.1141 
-0.3660 

\ -2.3243 



12.1560 
4.0252 

-1.1032 
4.5281 
0.0816 

-0.5332 

32.6228 
21.7486 
1.0356 
1.5535 
-0.2589 
-0.3884 



22.2650 
-11.7987 
19.2386 
42.7264 
-1.2588 
-3.5577 

-6.4667 
8.0604 
-9.6435 
-11.9660 
1.0987 
1.1326 

49.8089 
-16.6030 

1.8448 
-5.5343 
-0.4612 

1.3836 






-24.6846 
-11.7014 
2.7564 
2.9357 




10.6219 
-27.0825 
-1.0803 

2.7171 




-0.6250 
0.2246 
0.0887 
0.1026 






-12.9233 
-35.4784 
0.6168 
2.4965 




14.5052 
6.1964 
-1.6385 
-0.2564 






-4.0085 
-14.1041 
0.2854 
1.5568 




3.3472 
23.6695 
-0.2612 
-1.9148 



\ 


2.0820 

4.9828 
-0.2620 
-0.4249 j 






1.3651 
-4.8514 
-0.1847 
0.3886 



\ 






-6.6619 
1.9350 
1.1155 

-1.8214 






2.8566 
17.5058 
-0.1290 
-1.3503 






0.7622 
-0.7181 
-0.0387 

0.0524 



\ 
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/ o 











































4.8111 


-4.4336 


1.6880 


0.7207 








-12.2667 


-1.8940 


11.9366 


-2.5613 








-0.4893 


0.5008 


-0.1317 


-0.0975 


V 





1.2307 


0.0784 


—0.9657 


0.2051 


) 


/o 














\ 






























-1.3169 


-0.7444 


0.4827 


-1.2075 










0.4733 


0.2162 


2.9582 


1.1376 










0.1869 


0.1247 


-0.0218 


0.0613 




\o 





0.2161 


-0.2035 


-0.2282 


-0.0830 j 





(49) 



(50) 



/ - 



V 



1 

3 
_2 
9 
2 
189 
_J_ 
63 
1 
378 

1 
252 



_2 

3 

2 
9 

__2_ 

81 
_2_ 
27 

1 
162 
_ J_ 

54 






0.0129 
-0.0046 
-0.0018 
-0.0021 






0.0497 
-0.0144 
-0.0083 

0.0136 






-0.0092 
-0.0562 
0.0004 
0.0043 



/ 









\ o o 






0.7557 
-0.2716 
-0.1072 
-0.1240 






-0.1643 
0.0477 
0.0275 

-0.0449 






0.0861 
0.5277 
-0.0039 
-0.0407 






0.3224 
-0.3038 
-0.0164 

0.0222 



\ 



-0.0182 

0.0171 

0.0009 
-0.0012 J 

\ 



(51) 



(52) 



As far as the LO and NLO corrections parameterized by and are concerned 

our results agree perfectly with the findings of [11]. Contrariwise, the resummation of the 
NNLO logarithms is entirely new, and the corresponding matrices cf \ cf \ ef \ e^, 
and have never been computed before. 



7 Renormalization Scheme Dependences 

We would now like to elaborate on the question of renormalization scheme dependences in 
explicit terms, to gain an insight on how they arise beyond the LO, how various quantities 
transform under a change of scheme, and how the scheme dependences cancels out in 
physical observables. In this respect we will not only extend the existing NLO results 
[9, 10] to the NNLO level, but will also discuss the conceptual features related to the 
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renormalization of a s that, to the best of our knowledge, have not been studied in the 
context of the renormalization of effective field theories so far. 

It is well-known that beyond LO various quantities such as the ADM or the Wilson 
coefficients depend on the scheme adopted for the renormalization of the operators present 
in the effective theory. This scheme dependence arises because the requirement that all 
UV divergences are removed by a suitable renormalization of parameters, fields as well as 
operators, does not fix the finite parts of the associated renormalization constants. Indeed, 
these constants can be defined in different ways corresponding to distinct renormalization 
schemes, which are always related by a finite renormalization. In the framework of di- 
mensional regularization one example of how such a scheme dependence may occur is the 
treatment of 75 in n = 4 — 2e dimensions. In this context two well-known choices of scheme 
are the naive dimensional regularization scheme [33] with 75 taken to be fully anticom- 
muting and the 't Hooft-Veltman (HV) scheme [34] which comprises a 75 that does not 
have simple commutation properties with respect to the other Dirac matrices. Another 
example is the scheme dependence related to the exact form of the physical and evanescent 
operators chosen in the effective theory. We will discuss the latter issue in great detail in 
the following section. 

In order to show that physical quantities do not depend on the renormalization scheme 
and on the choice of the operator basis, we have to demonstrate how this dependence 
cancels out in the matrix elements of the effective hamiltonian introduced in Eq. (1) with 
C(n) given by Eq. (15). First, let us denote by 7^, and 7?, rjp with i — 1,2 the results 
for the ADM and the local parts of the matrix elements, that is, the finite pieces without 
logarithms of external momenta divided by the renormalization scale squared, obtained 
in two different renormalization schemes — see Eqs. (5) and (13). Furthermore, let us 
assume without loss of generality that the first scheme, which we shall call reference scheme 
hereafter, is distinguished from the other ones by the subsidiary condition rjp = = 0. 

It should be clear that for any given scheme a we can always switch to the reference 
scheme by the following finite renormalization: 



The corresponding transformations of the 0(af) and O(a^) anomalous dimensions is easily 
obtained using Eq. (31). At the NLO we reproduce the well-known result [9,10] 




(53) 



76' 



(i) 



7 «-[f« 7 (0) ]-2A)r, 



(i) 



(54) 



a 



whereas at the NNLO we find 



7* 2) =7?>" [r?>,7 (0) ] " K 1} ,7«] +r« [r«7 (0) ] - 4ftr?> - 2(3^ + 2 /? (f«) 2 . 



(55) 
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Similarly, for the transformation of the Wilson coefficients through NNLO we obtain 

T 



cm 



1 + 



4tt a + 



An 



r 



(2) 



(56) 



where again the NLO result is known for quite some time [35]. The relations connecting 
the ADM and the Wilson coefficients in two different schemes a and b can easily be derived 
from the above equations. 

With Eqs. (54) and (55) at hand, it is now straightforward to show that the matrices 
and introduced in Eq. (16) are independent of the renormalization scheme and the 
form of the operators considered. We start from the anomalous dimensions in the reference 



scheme and 7^ . These matrices can be accessed from any arbitrary scheme a using 
Eqs. (54) and (55). Let us transpose the latter equations and eliminate % and % by 
means of Eq. (A. 2). Finally, dropping the unnecessary subscript a, we obtain 



:,(2) 



Po 
Po 



fo)T £(i)l -2f3 R^, 



7 



(0)T #2) 



7 V 



(o)T #1)] + U(o)r^(i)l #1) 



(57) 



4/5 ^ (2) -2/5 lJ R«+2/5 (i? (1) ) 2 , 



which proves the scheme independence of M 1 ^ and R( 2 \ 

Next, A(°\ and A^ 2 \ obtained from the calculation in the full theory, clearly do not 
depend on the particular choice adopted for the renormalization of operators. In conse- 
quence, the factor to the right of LA°)(/i, /i ) in C(fJ,), as given in Eq. (15), which is related 
to the upper end of the evolution, is independent of the renormalization scheme. The same 
is true for the LO evolution matrix (/x, /x ). However, C(fx) still depends on the renor- 
malization scheme through K(fJ,) and consequently on and J^ 2 \ entering the Wilson 
coefficients to the left of U^(fj,, no). As is evident from Eqs. (7) and (13), this dependence 
on the lower end of the evolution is canceled by the one of the matrix elements (Q T (fi)). 
We have therefore explicitly seen that the matrix elements of the effective hamiltonian and 
the resulting physical amplitudes are scheme independent. 

It is important to emphasize that the renormalization scheme dependence discussed 
above refers to the renormalization of operators only, and has to be distinguished from 
the renormalization scheme dependence related to a redefinition of the strong coupling 
constant. In the following we will discuss the latter issue in detail, illustrating which 
effect a change in the charge renormalization has on the miscellaneous ingredients of the 
renormalization group improved perturbation theory. Finally, we will also prove that the 
matrix elements in the effective theory are invariant under a change of coupling constant. 
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The coupling parameter a s (p) in the new scheme will be denoted by a' s (p) and the 
finite renormalization relating the two schemes is written as 



a s (p) = 



a'(p) ( a'(p) . 



<*M • (68) 



It is now easy to show that the first two terms in the strong coupling expansion of the 
QCD j3 function, that is, Pq and fli, are scheme independent, while the third coefficient, 
namely /?2, transforms non-trivially [36]: 

& = + ci/3i - (ca - cj) p ■ (59) 

The renormalization constant matrices Z and Z' of the two different schemes can in 
general be related through a finite renormalization p(g') in such a way that 

°° / a' 2 \ k 

Z = p(g>)Z>, with pG/) = l + £ P (k) - (60) 

k=i ^ 77 ' 

Inserting this into Eq. (31) we obtain for the ADM in the primed scheme 

i'(g') = PigT 1 (i(g) + ^) M) ■ (6i) 

Comparing Eqs. (31), (53), (60) and (61) it should be clear that the NLO counterpart of 
Eq. (54) corresponding to the latter transformation reads 

y« = f i) _ [p«, f o)] _ 2/3o/) (i) + Ci ^ ( o) j (62) 



whereas the NNLO analog of Eq. (55) takes the following form 

f(2 ) = .(2) _ [p(2)^(o)j _ -4/3 p( 2 ) 
- 2/3 1 p( 1 ) + 2/3 (p^)) 2 - Cl [p« 7W] + c 27 {0) + 2^) . 



(63) 



Needless to say, that the finite operator renormalization p is unambiguously determined 
by the exact form of the change of charge renormalization parameterized by c\ and ci- 
Yet, it is not difficult to understand that by making use of the prior established fact that 
the matrix elements of the effective hamiltonian introduced in Eq. (1) are renormalization 
scheme independent one can always remove the dependence on and p^ from Eqs. (62) 
and (63). In other words, it is always possible to choose a scheme in which the operator 
renormalization constants remain pure UV poles. In consequence, the explicit relations 
between the coefficients codifying the finite charge and operator renormalization are not 
needed to prove the invariance of the effective theory under a redefinition of the strong 
coupling constant. Therefore we will not give this relations here. 
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It is easy to see, that the suitable change of renormalization scheme that one has to 
perform in order to remove the finite renormalization matrices p^ and p^ in Eqs. (62) 
and (63), is characterized through 

r« = -p«, and = -p^ + (p«) 2 , (64) 

and implemented by Eqs. (54) and (55). While the strong coupling constant is obviously 
invariant under such a change of renormalization scheme, that is, 

a':(p)=a> s (p), (65) 

the anomalous dimensions beyond LO transform non-trivially. Combining Eq. (60) with 
the latter transformation we obtain the following NLO anomalous dimensions in the double 
primed scheme 

y/(D = ^(i) + Ci ^(o) ? (66) 

whereas the NNLO result reads 

f{2) = ^(2) + c ^(0) + 2c ^(l) ) ( g 7) 

and the dependence on p^ and p^ has dropped out from the last two equations as in- 
tended. 

Next, let us write down the explicit expressions for the matrices J*- 1 -* and in the 
double primed scheme. Using Eq. (A. 2) and taking into account that the expansion co- 
efficients of the QCD (3 function in the double primed scheme coincide with those in the 
primed one, we find up to the NNLO level 



T"(l) _ 7(1) _ _£La(o)t 

2A, 7 ' 

J"(2) = j(2) + Ci J(l) _ iLj(^(0)T + ±Z121^)T + 



(68) 



with J« and defined as in Eq. (10). 

Furthermore, in order to separate the coupling scheme from the renormalization scheme 
dependence, let us assume without loss of generality, that the local parts of the matrix 
elements in the unprimed scheme fulfill the constraint f ^ = = 0. From Eqs. (66) and 
(67) it is then immediately clear that the local parts of the matrix elements in the double 
primed scheme correspondingly satisfy the relation r"^ = f"^ = 0. These two subsidiary 
conditions together with Eq. (68) now imply that the matrices and w 2 ^ in the double 
primed scheme are given by 



2A)" 7 



.2 .... :i . (69) 



2po 4p 8p 
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with R.W and denned as in Eq. (16). 

For the sake of completeness let us remark that only the third coefficient in the ex- 
pansion of the amplitude calculated in the full theory transforms non-trivially under a 
redefinition of the strong coupling constant: 

To prove that the effective theory does in fact not depend on the renormalization scheme 
employed for the strong coupling constant, let us first write down the analog of Eq. (15) in 
the double primed scheme: 



An 



X"W _ #"(i)^4"(o) 



+ 



An 



>(2) _ ^(1)^(1) _ (>(2) _ (£"(1)) 2 ) /"(O) 



(71) 



where 



and 



(72) 
(73) 

with V and Oj defined as in Eq. (9). 

It is now a matter of simply algebra to show that the LO evolution matrix in the double 
primed scheme is related to the LO evolution matrix in the unprimed scheme by 

f>"(°Wo) = L(/i)^ (0 Wo)£-Vo) , (74) 

where 



An 2[3 



An 



J _ ggOfo) c l >(Q)T , / gaQfo) 



4/3 7 8/^ 7 ' 

3C 1 ~ 2C2 J.(0) T C^(0) T\ 2 

4/3 7 8/^ 7 V 



(75) 



— # 

Taking into account Eq. (70) and expanding the double primed Wilson coefficient C"(fi) in 
terms of the unprimed coupling parameter leads to 



Os(/fo) 
47T 



Id) _ r»WaM 



+ 



An 



>) _ (id) _ Cl i(°)) - (>( 2 ) - (i^ 1 )) 2 ) a<°>" 



(76) 
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Inserting Eqs. (68), (69), (72) and (75) into the latter equation we finally obtain 
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+ 



An 



(77) 



>) _ ^(Dld) _ ^(2) _ (£(D) 2 ) 1(0) j > 



which shows that C"([/,) is nothing but C ((/,). Since does not depend on the scheme 
used to renormalize the high scale coupling constant, the same is obviously true for C"(fi). 
However, C"(/i) still depends on the renormalization scheme through K(fi) and conse- 
quently on and J^ 2 \ entering the Wilson coefficients to the left of LA°)(/x, // ). As is 
evident from Eqs. (7) and (13), and has already been discussed before, this dependence on 
the lower end of the evolution is canceled by the one of the matrix elements (Q T (fi)). We 
have therefore explicitly seen that the matrix elements of the effective hamiltonian and the 
resulting physical amplitudes are in fact invariant under a change of coupling constant. 



8 Change of Operator Basis 

In n = 4 dimensions, a change of the physical operators is always equivalent to a simple 
linear transformation 

Q' = RQ, (78) 

parameterized by a rotation matrix R, which affects the renormalization constants and the 
ADM in a trivial way: 

Z' = RZR-\ and 7' = R^RT 1 . (79) 

In the framework of dimensional regularization, the transformation corresponding to 
the change of basis turns out to be more complicated, as it generally involves evanescent 
operators as well. This feature basically reflects the fact that in order to formulate consis- 
tently the dimensional regularization of a theory containing fermionic degrees of freedom, 
the Dirac algebra has to be infinite-dimensional, which implies that evanescent operators 
are necessary to form a complete basis in n = 4 — 2e dimensions. In consequence, specify- 
ing the evanescent operators is necessary to make precise the definition of the MS scheme 
in the effective theory beyond LO, as can been seen for instance in Eq. (35). Clearly, 
EOM-vanishing operators are irrelevant to the present discussion. 

As long as the change of basis does not mix physical and evanescent operators, the 
ADM still changes in a trivial way. In particular, a linear transformation of evanescent 
operators does not affect the physical ADM at all. However, when the change of basis 
involves linear combinations of evanescent and physical operators, the situation turns out 
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to be more complicated [11]. Indeed, as we will explain in a moment, the new ADM is still 
given by Eq. (79), but the presence of evanescent operators induces a finite renormalization 
constant for the physical operators in the new basis. In order to restore the standard MS 
scheme definitions, a change of scheme is therefore required. 

Let us first consider a change of basis that consists of adding some evanescent operators 
to the physical ones, 

Q' = Q + WE, (80) 

parameterized by the matrix W. In this case the new ADM is still given by Eq. (79) 
because of the absence of mixing of evanescent into physical operators in the original basis. 
However, after the above transformation, the renormalization matrix corresponding to the 
physical operators in the new basis will contain a finite, non-vanishing contribution 

= WZ<$ , (81) 

where the subscript Q and E denotes an element of the physical and evanescent operators, 
respectively. In order to re-impose the standard MS conditions, the latter contribution 
must be removed by a change of scheme, implemented by Eq. (53). 

The situation is very similar for a change of basis that consists of adding multiples of 
e times physical operators to the evanescent ones 

E> = E + eUQ, (82) 

parameterized by the matrix U. In this case the ADM is unchanged because of its finiteness. 
However, the renormalization matrix of the physical operators in the new basis will contain 
a finite, non-vanishing contribution as well: 

zSf = -Z^U. (83) 

Needless to say, the above contribution must again be removed by a suitable change of 
scheme, in order to abide by the standard MS renormalization conditions. 

We therefore conclude in full generality that a change of basis in dimensional regulariza- 
tion is equivalent to a rotation plus a change of scheme. If we discount possible /i-dependent 
rotations of the operator basis, it should be clear from the discussion above that the most 
general change of basis comprises the three linear transformations of Eqs. (78), (80), and 
(82), as well as a rotation of the evanescent operators, which will be parameterized by the 
matrix M in what follows. In total we thus have 

Q> = r(q + WE} , and E' = M (eUQ + 1+e UW Ej . (84) 

The corresponding residual finite renormalization can be derived with simple algebra. Up 
to second order in a, we find 



7'(i.o) _ d 
QQ - K 



R- 1 
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J QQ 
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wz^ - (zg? + wz%§ - - \z^z%§w 



(85) 



With these expressions at hand, it is now straightforward to deduce how the ADM and 
the initial conditions of the Wilson coefficients transforms under the change of basis as 
given in Eq. (84). Up to the NNLO we obtain 
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(87) 



For what concerns the NLO parts of Eqs. (85), (86) and (87) our findings resemble the 
formulas derived in [11], if one takes into account that our definition of Z'qq^ differs slightly 
from the residual finite renormalization matrix used in the latter article. On the other hand, 
the complete NNLO relations Eqs. (85), (86) and (87) have never been presented before. 

After these general considerations, let us discuss in some detail how the anomalous 
dimensions given in Eqs. (36), (37) and (38) are transformed in going to the "traditional" 
basis of physical operators [1, 10, 14] 



Qi 

Qa 



(*l7^)E,(^&, 



(88) 



In the above definitions a and (3 denote color indices. 

The one- and two-loop evanescent operators that accompany the "traditional" basis 
can be found by imposing the requirements given in [14]. At the one-loop level they are 



E'i 



'(i) 



(^Ll m ciMl m bl) ~ (16 - 46) Q[ , 
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E' {2 
Ef 

eT 



E' 2 {1) = (sll^ 2 ^)(4l m 4) - (16 - 46) Q' 2 , 
Ef ] = {sl lm bl) J2 " (16 - 46) Q' 3 , 

EP = (s a L ^ m bl) J2 q (& m <l a L) - (16 - 46) Qi , 
E' 5 {1) = {sl lm bt) E ? (^™^) - ( 4 + 4e ) > 

= (s a L l m bl) £ t^™^) - (4 + 46) . (89) 
Following the same procedure, we find the following two-loop evanescent operators: 

E'i 2) = (^7, 1 M 2 .3M4, 5 c?)(c?7 Ml ' i2M3M4At5 ^) " (256 - 224 e ) Q[ , 
= (^7, 1 M 2 M3M4M 5 c2)(c?7 w ^' 13 ' i4M5 &?) - (256 - 2246) Q' 2 , 
= (sll m ^M) Y, q ^ m ^ 5 4) ~ (256 - 2246) Q' 3 , 
= {sll m ^A) Y. q ^Ll mm Ql) - (256 - 2246) Q' A , (90) 
= £^7^™^) " (16 + 1286) Q> 5 , 

{sll m ^A) £^7™"^) - (16 + 1286) Q> 6 . 

It turns out that in order to transform the ADM given in Eqs. (36), (37) and (38) from 
the initial set of operators to the "traditional" basis, we have to introduce four additional 
one-loop evanescent operators 

eP = (sliM E g (^ M ^) - jjCs + \qb , 

Ejp = (s L ^T a b L ) £ (qi^T a q) - \q a + \q, , 

^— '9 3 

£? } = {sLi m b L ) £ g (gV™ 3 7 5 g) - f g 3 + \q, , 

as well as four additional two-loop evanescent operators 

Ef ] = {s Llm ^b L ) J2 q (qr m ^ 5 l 5 q) ~ ™Qs + f Q 5 , 
Ef ] = (s L l^ 2m T a b L ) J2 q (qi^ m5 l 5 T a q) - ™Q 4 + ^Q 6 , 
Ef = (s Ll , m4m b L ) J2 q (qi mm ^l 5 q) ~ ^Qs + , 
E% ] = {s Llmm , 7 T% L ) J2 q (qi mm ^T a l5 q) - ^Q 4 + ^Q 6 . 
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(91) 



(92) 



It should be clear, that the evanescent operators E^-Ejp and Ejp-Ejp are not needed as 
counterterms in the initial basis of operators. However, some linear combinations of them 
will become parts of either the physical or the evanescent operators in the "traditional" 
basis through the change of basis given by Eq. (84). 

At this point a comment concerning the computation of the renormalization constants 
involving the insertions of the additional evanescent operators is in order. Transforming 
the three-loop anomalous dimensions from the initial to the "traditional" basis requires the 
knowledge of one- and two-loop diagrams with insertions of E^—E^p, which introduces 
traces with 75 into the calculation. In this context we follow [37], and avoid anticommu- 
tation of 75 in any fermionic line containing an odd number of 75. Moreover, while traces 
containing an odd number of Dirac matrices and a single 75 do not pose any problem as 
they vanish algebraically, we do not evaluate traces containing an even number of Dirac 
matrices and a single 75 in n = 4 — 2e dimensions. This brings to life new evanescent 
operators — we will call them trace evanescent in the following — that can in general be 
written as a contraction of a suitable Dirac structure with one of the following evanescent 
tensors 

= Tr (7mi...m™7s) - Tr (7 w ...^ m 75) , (93) 

where the four-dimensional traces Tr (7 Atl ... Mm 7s) can be calculated recursively from the 
initial value Tr (75) = applying 

m 

(7/ui.../im75) = ^ ] ( — 1) g^ 1>lj Tr (7// 2 — ^-1^+1 — ^75) — g 6 ^ 1 1 2 3 Tr (y^... ^ mUl u 2 u 3 ) ■ (94) 
i=2 

Here g^ 2 = diag (1, —1, —1, —1) denotes the four-dimensional metric tensor, and e^ lfM2fM3 ^ 4 
is the totally antisymmetric Levi-Civita tensor defined so that e 0123 = 1. Furthermore, the 
second trace in the above equation should also be taken in n = 4 dimensions. It can thus 
be computed recursively from Tr(l) =4 using 

m 

ft(7,i...O = 5Z(-l) i sr w ^(7 M ... w _i W+ i...p ra ) • (95) 
i=2 

Apparently, the trace evanescent operators originating from Eq. (93) have to be treated 
on the same footing as the regular ones introduced earlier on. The idea of introducing 
more and more evanescent operators seems to make the use of an naive anticommuting 75 
in multi-loop calculations involving chiral operators futile. Fortunately, for the problem 
at hand this is not the case, as it turns out that the one-loop insertions of E^ and 
E^ needed to find the transformation of the two-loop anomalous dimensions between the 
initial and the "traditional" basis involve only the trace Tr(7 MlAt2 75), which however is zero. 
This observation has been made already in [11]. Furthermore, to find the transformation 
of the three-loop anomalous dimensions matrices, both one-loop insertions of E^-E^ 
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and two-loop insertions of E\ and Eq are required. Typical examples of non-vanishing 
two-loop 1PI diagrams are shown in Figure 6. However, also in this case the number of 
new evanescent structures is rather small, since the necessary operator insertions introduce 
only the non-trivial trace Tr(7 MlM2M3M4 7 5 ). The complete list of trace evanescent operators 
relevant to find the transformation of the three-loop anomalous dimensions between the 
initial and the "traditional" basis is given in Appendix A. 2. 

If possible any regularization prescription should respect all symmetries of the bare 
theory, such as gauge and BRST invariance encoded in the Ward and Slavnov- Taylor 
identities, or Bose symmetry. It is interesting to note that all this requirements are not 
necessarily fulfilled for an arbitrary choice of trace evanescent tensors. For example, adding 
any multiple of e times the four- dimensional traces Tr (7 Ml ... Mm 75) to the right hand side of 
Eq. (93), in general, spoils the usual Ward and Slavnov- Taylor identities as well as the Bose 
symmetry of some of the resulting 1PI Green's functions even after the correct subtraction 
of all subdivergences. In order to verify that neither problem arises with the definition of 
trace evanescent operators adopted in Eqs. (A. 4) and (A. 5), we have calculated the pole and 
finite parts of the off-shell b — > sqq, b — * s, b — > sg, b — * sgg and b — > sggg matrix elements 
with one-loop insertions of E^p-E^ and two- loop insertions of Ejp-EP, and checked 
explicitly that i) the resulting operator renormalization constants are independent of the 
external states used in the calculation, and that ii) the subtracted 1PI Green's functions 
with two and three external gluons are cyclic under the interchange of any two gluons as 
required by Bose symmetry. Concerning the former issue, let us mention, that in order to 
decompose the finite parts of the two-loop matrix elements of eP~eP corresponding to 
the b — > sgg and b — > sggg transitions, one has to enlarge the off-shell operator basis to 
contain besides the EOM-vanishing operators nP and N^-N$ the following four gauge 
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non-invariant operators of dimension-six: 

A? = e mtli G a ^d^G a ^{s Ll ^b L ) , 

4 2) = e, m J abc G b ^d^G^\rs Ll ^T a b L ) , 

ro\ (96) 
4 2) = e mz ,J abc G a ^G b ^G c ^{s Ll ^b L ) , 

Af = e MTOM4 (d abe f cde + d ace f dbe + d ade f bce ) G a ^G b ^G c ^{s Ll ^T a b L ) . 

Of course the finite two-loop renormalization between E^-E^ and the anomalous two- 
and three-gluon operators A^-A^ does not affect the residual finite renormalization ma- 
trices Z'qq^ and Z'qq^ defined in Eq. (85) which are unambiguously fixed by the corre- 
sponding one- and two-loop on-shell b — > sqq and b — > sg matrix elements. In consequence, 
the presence of the gauge non-invariant operators A^-A^ leaves the final result for the 
anomalous dimensions of the four-quark operators Q\-Qq in the "traditional" scheme un- 
altered up to the three-loop level. 

In case the above considerations might not fully convince a suspicious reader that our 
treatment of traces containing an odd number of Dirac matrices and a single 75 is con- 
sistent, it may be worthwhile to justify it in another independent way. To exclude all 
possibility of doubt concerning our regularization scheme, we have computed the renor- 
malization constants involving the insertions of Etp-E^ using the HV definition of 75 in 
n = 4 — 2e dimensions, and verified explicitly that, for what concerns the non-anomalous 
operators, the latter results coincide with those obtained by the "dimensional reduction" - 
like treatment of traces containing 75, implemented by Eqs. (93), (94) and (95). For the 
sake of completeness, let us also mention, that in order to decompose the finite parts of the 
two-loop b — > sgg and b — > sggg matrix elements with insertions of £5 and Eq computed 
in the HV scheme, not four but eight anomalous operators are required. Of course, the ap- 
pearance of these additional gauge non-invariant operators, which correspond to Af^-A^ 
with right instead of left chiral quark fields, does not alter the final result for the three-loop 
0(al) anomalous dimensions in the "traditional" scheme. 

As regards the relevance of the trace evanescent operators defined in Eqs. (A. 4) and 
(A. 5), let us finally point out, that it would be incorrect to conclude that they are only 
needed to find the transformation of the three-loop anomalous dimensions matrices be- 
tween the initial and the "traditional" basis of operators. In fact, their role in the "tra- 
ditional" basis is the same as the one of the evanescent operators introduced in Eq. (90), 
as they, together with the latter ones, make precise the definition of the renormalization 
scheme in that basis at the NNLO. As should be clear from what has been said in the last 
section on the cancellation of scheme dependences in general, the corresponding scheme 
dependence is canceled by the one of the two-loop 0(a 2 s ) matrix elements, which have to 
be calculated using exactly the same definition of trace evanescent operators. 

The renormalization constant matrices entering Eq. (85) are found from one- and two- 
loop matrix elements of physical and evanescent operators. We give the relevant ones, as 
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well as the matrices characterizing the change of basis in Appendix A. 3. Our final results 
for the residual finite renormalization read 
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and 
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While the one-loop residual finite renormalization Z'qq^ is independent on the choice of 
trace evanescent operators, the / parts of the two-loop residual finite renormalization 



Z'qq" 1 relating the QCD penguin operators in the unprimed and primed scheme are not. 



The / parts given in Eq. (98) correspond to the specific choice of trace evanescent operators 
adopted in Eqs. (A. 4) and (A. 5). 

With these expressions at hand, it is now a matter of simple algebra to find the ADM 
and the initial conditions of the Wilson coefficients in the "traditional" basis. Using 
Eqs. (36), (37), (38), (86), (97) and (98) we obtain for the regularization- and renormaliza- 
tion scheme independent one-loop 0(a s ) anomalous dimensions matrix: 
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-16+f/ / 

While the matrix 7^°) is renormalization-scheme independent, 7^ 1 - ) and 7'^ are not. In 
the MS scheme supplemented by the definition of evanescent operators given in Eqs. (89) 
and (90) we obtain for the two- loop 0{a 2 s ) anomalous dimensions matrix: 
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The / and f 2 parts of the matrix -y^ 2 ) describing the mixing of the QCD penguin operators 
do in addition depend on the definition of trace evanescent operators given in Eqs. (A. 4) and 
(A. 5). In the corresponding scheme we find for the three-loop O(af) anomalous dimensions 
matrix: 



f( 2 ) = 



19859 1543 f , 298 f 2 , 80 r /- 
36 18 J ' 81 J ' 3 J ^ 3 



4741 
12 



f /-§/ 2 -(672+80/)C 3 








f^f/-^/ 2 -(672+80/)C3 

1165 7 f , 82 f 2 , 80 f r 
36 18 J ~ + ~ 81 ■> ' 3 J ^ 3 








52488 
22475861 , 1025695 

13122 t 17496 
146039 , 8338543 r_ 
81 ' 52488 J 
100832 



58231 | 595 r 
1944 + 81 ^ 
3664721 , 23831 f , 280 f 
' < 4374 J * 81 S3 



/+S/ 2 

15961 



560 
81 



f/)C3 
/ 2_( 672+ 6200 / ) C3 



17705 r J_ r2 
J 27 J 



81 486 J 27 

10779689 t i 18005 i-2 , 280 



> f I 18005 r2 
~J ' 4374 J 



17496 
12604483 1382815 



77239 | 53 r 
648 '27 J 
23111^4024^ 



1458 
7 



27 

26192 



+80/) Cs 



31679 2583223 f 6359 t 2 3304 j ^ 



27 17496 ■' 1458" 

60448 i 25841 r 



81 

128992 , 7661297 f 31109 f 2 5464 
81 " r 17496 J 1458 J 27 



81 
23 /2 



9 



/Cs 



(101) 



261287 , 163 
1944 ' 81 
2612539 i 38285 



f , 20 a 
j I 7228 

J ~T 01 



13122 1 17496 
437203 1360099 



1049149 f _49 f 2 , ( : 
17496 ■< 27' """V" 



324 
1826987 



r , 13289 r2_ 

26244 ' + 4374 J 



V)C3 

;f+ 2 i i /)c3 



T348 972- i /+i/ 2 -(378+20/)C3 

1091127 , 509723 t , 27815 c 2 , ^210 I 7228 /)^3 



1944 



26244 



f I 27815 f 2 
J ' 4374 •< 



648 1 27 
1722187 16541 



+ P/-20C3 



4374 5832 J ' 9 J 



8748 
266077 



L C3 



199123 i 162707! 



108 1 8748 
472667 . 88933 
216 ' 324 
1946849 , 5194621 



.(4088 +2Q/ ) C3 

^-fif/ 2 -(40+2i4/)C 3 
/-li/ 2 +(462+60/)C3 
/-f!i/ 2 +(378 + 2 #/)C3 J 



As far as the one- and two- loop self- mixing of the four-quark operators Q[-Q'q, namely 
7'(°) and 7'^ are concerned, our findings agree perfectly with the results of the direct 
computations [10,14]. Furthermore, they also agree with the results presented in [11], 
which were obtained by performing a change of scheme from the unprimed set of operators 
to the primed one. On the other hand, the three-loop self-mixing of Q' x -Q' % described by 
7'( 2 ), is entirely new. 

Similarly, employing Eqs. (26), (87), (97) and (98) we find for the initial conditions of 
the Wilson coefficients in the "traditional" basis up to 0(a 2 ): 
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Whereas the two-loop 0(0;^) corrections to the initial conditions C[(M W )-C'§(M W ) are 
entirely new, our findings for the one-loop 0(a s ) corrections agree again perfectly with all 
preceding calculations [10,11,14]. Finally, let us mention that the 0(a 2 s ) corrections to the 
initial conditions C[(M W )-C' 6 (M W ) do not depend on the special choice of trace evanescent 
operators made in Eqs. (A. 4) and (A. 5). The NNLO analytic formulas for the low-energy 
Wilson coefficients relevant for non-leptonic B meson decays in the "traditional" scheme 
can be found in Appendix A. 4. 

9 Summary 

In this paper we have extended the SM analysis of the effective hamiltonian for non-leptonic 
|AF| = 1 decays to the NNLO. The main ingredient of this generalization is the three-loop 
ADM describing the mixing of the current-current and QCD penguin operators, which we 
have computed in an operator basis that allows to consistently use fully anticommuting 75 
in dimensional regularization to all orders in perturbation theory. The issue of renormaliza- 
tion scheme dependences, their cancellation in physical quantities, and the transformation 
properties of the ADM and the initial conditions of the Wilson coefficients under a change 
of scheme has been discussed thoroughly. In particular, we have elaborated on the scheme 
dependence related to the renormalization of the strong coupling constant, a feature that, 
to the best of our knowledge, has not been studied in the context of the renormalization 
of effective field theories so far. As a practical application of our general considerations, 
we have derived the explicit NNLO relation between our and the so-called "traditional" 
renormalization scheme, which allowed us to calculate indirectly the NNLO ADM and 
the corresponding matching conditions in the latter scheme. Finally, we have solved the 
RGE to obtain the analytic expressions for the low-energy Wilson coefficients relevant for 
non-leptonic B meson decays through NNLO in both schemes. 
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Appendix 

A.l Derivation of and 5< 2 ) 

In order to derive the explicit expressions for the matrix kernels S^ 1 ) and as given 
in Eq. (11), we follow [9,39] and compute the partial derivative of Eqs. (4) and (6) with 
respect to g. After some algebra one finds the following differential equation for K(g) 
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Inserting Eq. (7) into the last equation we obtain 
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for the parts proportional to g and g 3 , respectively. After diagonalizing these equations 
with the help of Eq. (9) we find 
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Finally, solving the first equation for G^f and inserting the result into the second equation, 
one obtains the expression for the elements of as given in Eq. (11). 
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A. 2 Trace Evanescent Operators 

In the following we specify the exact form of the so-called trace evanescent operators arising 
from the one- and two-loop diagrams with insertions of E^p-E^p. At the one-loop level 
the specific structure of only one of them is needed: 

E? = {si mz P L b) Tr (y™^ 75 ) _ 24(s L7 ^6 L ) , (A.4) 
while at the two-loop level we encounter eight additional trace evanescent operators: 
E[ 2) = (s^P L b) ^(yr^g) Tr ( 7 ^<~ 7s) _ 4 0(f L7w 6 L ) J^ffy" 1 ?) 

+ ±{sLl m b L ) J2& m q) , 
E { ? = (si^PLl^b) Tr ( 7 ^~ 7s) _ 2 A{s Ll ^b L ) , 

Ei 2) = (s^P L ^ 3 b) Tr (y*^«7fe) + 8(^ 2 4 (s l7 ^l) - ^ 2 B (*l7^l)) , 

£? } = (s^PLj^b) Tr (y™/,4 75) + 8 ((s l7 ^ r ) - ^ 4 (*l&«)) , 

Ef ] = {s lm P L b) Tr ( 7 ™m B75 ) _ 8 (( S - L7 ^fe L ) - ^fa/frA)) , (A 5) 

= (s 7mi P l 6) Tr ( 7 ^3„4 75 ) + 4 ^ L7 M 2 ^3M4 6i) 

- g^ 2 ^ 3 (s L Y 4 b L ) + g^ 4 (s L Y 3 b L ) - ^ 4 (sj/T&l)) , 

= {s lm P L b) Tr (y^«7 B ) + %((s Ll ^b L ) - 2g£(8 L -f*b L ) 

+ 2~g£(s Ll ^b L ) - ~ 9 ^ 5 (sliM) , 
Ef ] = (s^P^^b) Tr + 8((s L %T 5 b L ) - 2~g£{s Ll ^b L ) 

+ 2~g£{s Ll ^b L ) - ~g^(SLlM) , 

where we have used the definitions T^.^j"^™ = 7/ m • • • 7^-7 w+1 • • • 7^™ an d -Pl = (1 — 7s)/2- 
Finally, let us recall two important features of the latter operators. First, their role in the 
"traditional" basis is the same as the one played by the evanescent operators in Eq. (90), 
that is, they are needed to define the renormalization scheme in that basis at the NNLO. 
Second, the / and f 2 parts of the three-loop ADM 7 '( 2 ) describing the mixing of the QCD 
penguin operators in the "traditional" basis depend on the definition of these operators. 
The above definitions correspond to the "dimensional reduction" -like treatment of traces 
with an even number of Dirac matrices and a single 75, employing Eqs. (93), (94) and (95). 
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A. 3 Change to the "Traditional" Operator Basis 



In order to give the explicit expressions for the matrices R, W, U and M characterizing 
the change to the "traditional" basis, we first have to define the unprimed and primed set 
of operators according to Eq. (84). The physical and evanescent operators in the initial 
basis are given by 



Q T — (Qi, ■ ■ ■ , Qq) , 
E T = {E^\ . . . , E^\ E[ 2 \ . . . , E^) , 
while the "traditional" basis consists of the following two sets of operators: 



(A.6) 



Q 

E' 



(Qi,- 



Qe 



p'W i/(2) 



p/(2) p (2) p (2) p (2) p(2)x 
^6 ?-^3 ?-^4 5-^7 5-^8 J 



(A.7) 



Needless to say, that E^-E^ and E^-E^ play the role of extra, in principle unnecessary 



operators in the initial operator basis. The same is true for Ef\ Ef\ E^ and E^' in 
the "traditional" basis. They are just included for completeness in the above equations. 
Although the trace evanescent operators introduced in Eqs. (A. 4) and (A. 5) are needed 
to make precise the definition of the "traditional" renormalization scheme we refrain to 
include them in Eq. (A.7), since they influence the transformation between the unprimed 
and primed set of operators only in an indirect way, that is, through the finite parts of 
the two-loop renormalization constants of some of the evanescent operators included in the 
initial basis. This will be explained at the very end of this appendix. 

With this definitions at hand, it is just a matter of simple algebra to find the explicit 
expressions for the matrices R, W, U and M . The rotation matrix R, which links the 
physical operators together, is given by 
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The matrix W parametrizes a redefinition of the physical operators Q by adding some 
evanescent operators E to them. In the case at hand, W reads 
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(A.9) 
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On the other hand, U describes a redefinition of the evanescent operators E by adding 
some multiples of e times physical operators Q to them. The relevant matrix U takes the 
following form 
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Finally, the matrix M represents a simple linear transformation of the evanescent operators. 
In our case we find 
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Parts of the above matrices have already been given explicitly in [11], where the change 
of basis from the initial to the "traditional" basis has been performed including the NLO 
corrections. If we take into account that the definition of E^-E^ adopted in Eq. (91) 
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differs slightly from the definition of E^-E^ used in [11], our results agree with the 
expressions given in the latter paper. 

The renormalization constant matrices entering Eq. (85) are found from one- and two- 
loop matrix elements of physical and evanescent operators. In the following we will give only 
the relevant entries of the necessary renormalization constant matrices, denoting elements 
that do not affect the final results for the residual finite renormalizations introduced in 
Eq. (85) with a star. For the finite renormalization between evanescent operators E and 

— * 

physical operators Q we get 
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For the one-loop mixing of physical operators Q into evanescent operators E we obtain 
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The one-loop mixing among evanescent operators E reads 
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At the two-loop order we find for the finite renormalization between evanescent operators 
E and physical operators Q: 
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— * — * 

The two-loop mixing of physical operators Q into evanescent operators E is given by 
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Finally, the two-loop mixing among evanescent operators E reads 
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As far as the one-loop renormalization constant matrices are concerned, let us note, that 
our results agree with the findings of [11], after taking into account that the definition of 



Ecp—E^ adopted in Eq. (91) differs slightly from definition of E^'-Eg 1 ' used in the latter 
article. On the other hand, the two-loop renormalization constant matrices involving the 
insertion of E^ and E^ , are entirely new and have to our knowledge never been computed 
before. Finally, let us mention that while the matrices Zq£\ Zqq\ Zq£\ Z^^ 

do not depend on the special choice of trace evanescent operators used throughout the 
calculation, the / parts of Z^n do depend on the latter choice. The / parts of Z^'q as 
given in Eq. (A. 15) correspond to the specific choice of trace evanescent operators adopted 
in Eqs. (A.4) and (A. 5). 



(i) 77.(1) 
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A. 4 Wilson Coefficients in the "Traditional" Operator Basis 

In this appendix we shall use the ADM given in Eqs. (99), (100) and (101) to find the 
explicit NNLO expressions for the Wilson coefficients in the "traditional" basis: 

CM = C,'<V) + >(*) + (^) V W , (A.18) 

with % = 1-6. Using the general solution of the RGE given in Eq. (6), we arrive at 
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(A.19) 
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(A.21) 



(A.22) 
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As far as the LO and NLO corrections parameterized by c v , c ( 



,'(o) 



(A.31) 



and are 



concerned our results agree perfectly with the findings of [9,10,35]. Contrariwise, the 
resummation of the NNLO logarithms is entirely new, and the corresponding matrices 



~/(2) ~/(2) ,7(2) ~/(2) ./(2) ;/(2) 



e' 1 lz; , e^' , t'^' and have never been computed before. Finally, let us 



/(2) 



7(2) 
c 2 ) 



?'(2) 



A 2 ) 
e 2 1 



i'^ and do not depend on the 



mention that while the matrices c 1 
special choice of trace evanescent operators used throughout the calculation, the entries of 
c describing the evolution of the Wilson coefficients of the QCD penguin operators do 
depend on the latter choice. The numbers given in Eq. (A. 25) correspond to the specific 
choice of trace evanescent operators adopted in Eqs. (A. 4) and (A. 5). Recalling what 
has been said earlier on the cancellation of scheme dependences in general, it should be 
clear, that the associated scheme dependence is canceled by the one of the two- loop 0(a 2 s ) 
matrix elements, which have to be calculated using the same definition of trace evanescent 
operators. 
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